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The ratio of the measured B —» D^lv decay rates for < = t vs. e,(i deviate from the standard 
model (SM) by about 4a. We show that the data are in tension with the SM, independent of 
form factor calculations, and we update the SM prediction for B(B —> X c rv)/B{B —¥ X c lv). We 
classify the operators that can accommodate the measured central values, as well as their UV 
completions. We identify models with leptoquark mediators that are minimally flavor violating in 
the quark sector, and are minimally flavor violating or r-aligned in the lepton sector. We explore 
experimental signatures of these scenarios, which are observable in the future at ATLAS/CMS, 
LHCb, or Belle II. 


I. INTRODUCTION 


Measurements of the B —> Dtv and B —► D* tv decay 
rates are now available from BaBar [2 d] and Belle 0 
with their full datasets. The B —> D*tv decay mode 
was also observed recently by LHCb [3J. These measure¬ 
ments are consistent with each other and with earlier 
results El IH, and together show a significant deviation 
from Standard Model (SM) predictions for the combina¬ 
tion of the ratios 


R(X) 


B(B ->• Xtv) 
B{B -± Xlv) ’ 


( 1 ) 


where l = e, /i. The measurements are consistent with 
e//i universality (7|, iSj. The R{D M) data, their aver¬ 
ages [55] , and the SM expectations are summarized 

in Tablc[I] (If the likelihood of the measurements is Gaus¬ 
sian, then the deviation from the SM is more than 4cr.) 
Kinematic distributions, namely the dilepton invariant 
mass q 2 , are also available from BaBar and Belle Eua, 
and must be accommodated by any model that modifies 
the rates. In the future, Belle II is expected to reduce 
the measured uncertainties of R(D^) by factors of ~ 5 
or more [T3] . thereby driving experimental and theory 
precision to comparable levels. 

In the type-II two-Higgs-doublet model (2HDM), the 
B -t tv rate (as well as B —» tv) receives contribu- 




R(D) 



R(D*) 


Corr. 

BaBar 

0.440 

± 0.058 ± 

0.042 

0.332 

± 0.024 ± 

0.018 

-0.27 

Belle 

0.375 

± 0.064 ± 

0.026 

0.293 

± 0.038 ± 

0.015 

-0.49 

LHCb 




0.336 

± 0.027 ± 

0.030 


Exp. average 

0.391 

-H 

o 

d 

-H 

0.028 

0.322 

-H 

GO 

o 

o 

-H 

0.012 

-0.29 

SM expectation 

0 

300 ± 0.010 

0 

252 ± 0.005 


Belle II, 50/ab 


±0.010 



±0.005 




TABLE I. Measurements of R(D ( *') [I] 3] [2, their aver¬ 
ages 0, the SM predictions [MS!, and future sensitiv¬ 
ity [13] • The first (second) experimental errors are statistical 
(systematic). 


tions linear and quadratic in mb m T tan 2 ft /m 2 H ± [HHE], 
which can be substantial if tan/3 is large. However, the 
R(D^) data are inconsistent with this scenario jTj. 

Discovering new physics (NP) in transitions between 
the third and second generation fermion fields has long 
been considered plausible, since the flavor constraints are 
weaker on four-fermion operators mediating such transi¬ 
tions. (Prior studies of B —> X s vv [T7j and H( s ) —> 
t + t~{X) [18[ 19] decays were motivated by this con¬ 
sideration.) However, B —i D^tv is mediated by the 
tree-level b —> c transition. It is suppressed in the SM 
neither by CKM angles (compared to other B decays) 
nor by loop factors, with only a modest phase space sup¬ 
pression due to the r mass. This goes against the usual 
lore that the first manifestations of new physics at low 
energies are most likely to occur in processes suppressed 
in the SM. 

The goal of this paper is to explore flavor structures 
for NP capable of accommodating the central values of 
the R(D^) data summarized in Table [i] To do so, a 
sizable NP contribution to semileptonic b —>■ c decays 
must be present, and the NP mass scale must be near 
the weak scale. This requires nontrivial consistency with 
other constraints, such as direct searches at the LHC and 
precision electroweak data from LEP. When NP cou¬ 
plings to other generations are present, constraints from 
flavor physics, such as meson mixing and rare decays, 
also play a role. For example, any flavor model predicts 
some relation between the be vt and bu vt operators, so 
models explaining R(D W) must accommodate the ob¬ 
served B~ —> tv branching ratio, which agrees with the 
SM 2IB |22. We show below that despite strong con¬ 
straints some scenarios remain viable and predict signals 
in upcoming experiments. 

We begin by presenting new inclusive calculations that 
demonstrate that the measured central values of R(D^) 
are in tension with the SM, independent of form factor 
computations. Then, in Sec. [TT| we perform a general 
operator analysis to identify which four-fermion opera¬ 
tors simultaneously fit R{D) and R(D*). In Sec. Ill we 


discuss possible mediators that can generate the viable 
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operators. We identify working models with leptoquark 
mediators that are minimally flavor violating in the quark 
sector, and we confirm their consistency with current ex¬ 
perimental constraints. Finally, Sec. |IV| contains our con¬ 
clusions and a discussion of possible future signals at the 
LHC and Belle II. Appendix |A| contains a discussion of 
U( 2) 3 models. 

A. Standard Model considerations 

The tension between the central values of the R(D 1*)) 
data and the SM is independent of the theoretical pre¬ 
dictions for quoted in Table 0 The measured 

R(D^) values imply a significant enhancement of the 
inclusive B —> X c tD rate, which can be calculated pre¬ 
cisely in the SM using an operator product expansion, 
with theoretical uncertainties that are small and essen¬ 
tially independent from those of the exclusive rates. 

To see this, note that the isospin-constrained fit for the 
branching ratios is quoted as jTl 

B(B -a D*tv) + B{B -a Dtp) = (2.78 ± 0.25)% , (2) 

which applies for B ± decays (recall the lifetime difference 
of B ± and B°). The averages in Table [I] imply for the 
same quantity the fully consistent result, 

B{B -a D*tD) + B{B -A Dtp) = (2.71 ± 0.18)% . (3) 

The SM prediction for R(X C ), the ratio for inclusive 
decay rates, can be computed in an operator product ex¬ 
pansion. Updating results in Refs. p21 [25] . and including 
the two-loop QCD correction [21], we find 

R(X C ) = 0.223 ±0.004. (4) 

The uncertainty mainly comes from m \ s , the HQET ma¬ 
trix element Ai, and assigning an uncertainty equal to 
half of the order a 2 term in the perturbation series in 
the IS 1 scheme 2B:. The most recent world average, 
B(B~ -A X c ev) = (10.92 ± 0.16)% [26]|27], then yields 
the SM prediction, 

B(B~ -A X c tp) = (2.42 ± 0.05)% . (5) 

In B~ -a X c ev decay, hadronic final states other than 
D and D* contribute about 3% to the 10.92% branching 
ratio quoted above, and the four lightest orbitally excited 
D meson states (often called collectively D**) account for 
about 1.7%. Using Ref. [21] for the theoretical descrip¬ 
tion of these decays, taking into account the phase space 
differences and varying the relevant Isgur-Wise functions, 
suggests R(D**) > 0.15 for the sum of these four states. 
This in turn implies for the sum of the central values of 
the rates to the six lightest charm meson states 

B{B->dWtD) + B(B->D**tv)~ 3%, (6) 

in nearly 3 a tension with the inclusive calculation in 
Eq. ©.' Note that Eqs. ([2|, Q, and ([6| are also in mild 


tension with the LEP average of the rate of an admixture 
of b -flavored hadrons to decay to r leptons [29], 

B(b -A Xt+p) = (2.41 ± 0.23)% . (7) 

Since both the experimental and theoretical uncertainties 
of B{B -a X c tp) are different from the exclusive rates, 
its direct measurement from Belle and BaBar data would 
be interesting and timely [30] . 


II. B -A D^tP operator analysis 

In this section we study operators mediating b -A ctP 
transitions. In contrast to prior operator fits ESHU, 
we adopt an overconrplete set of operators corresponding 
to all possible contractions of spinor indices and Lorentz 
structures to help with the classification of viable models. 
(We also take into account the constraints from q 2 spec¬ 
tra, which were unavailable at the time of the first oper¬ 
ator analyses.) Although Fierz identities allow different 
spinor contractions to be written as linear combinations 
of operators with one preferred spinor ordering, the set 
of possible currents that can generate the operators is 
manifest in the overcomplete basis. 

We parametrize the NP contributions by 

n = ^ u cfc o Vl + A E c i'" ] 0 ^ n) ■ ( g ) 

v i 

(Throughout this paper we do not display Hermitian 
conjugates added to interaction terms as appropriate.) 
Here the primes denote different ways of contracting the 
spinors, as shown in Table [Tl] which also presents their 
Fierz transformed equivalents in terms of the “canoni¬ 
cally” ordered fields (unprimed operators). In the SM, 
only the Oy r operator is present. (For illustration, 
the type-II 2HDM generates the operator Os R with 
C Sr / A 2 = — 2\/2 GpVcb mb m T tan 2 f3/m 2 H± .) 

We do not consider the possibility of the neutrino be¬ 
ing replaced by another neutral particle, such as a sterile 
neutrino, which yields additional operators. The large 
enhancement of an unsuppressed SM rate favors NP that 
can interfere with the SM. A non-SM field in the fi¬ 
nal state would preclude the possibility of interference, 
leading to larger Wilson coefficients and/or lower mass 
scales for the NP, making the interpretation in terms of 
concrete models more challenging. 

We assume that the effects of NP can be described 
by higher dimension operators respecting the SM gauge 
symmetries. This is only evaded if the NP mediating 
these transitions is light or if it is strongly coupled at 
the electroweak scale; in either case there are severe con¬ 
straints. We classify operators by the representations 
under SU(3)c x SU(2)l x U(1)y of the mediators that 
are integrated out to generate them, as shown in the last 
column of Table [H] Some mediators uniquely specify a 
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Operator 

Fierz identity 

Allowed Current 

<S£int 

Ov L 

Ov R 

Os R 

(c7 m Pl&) (f7 m Pl^) 
(c7 m P«6) ( t^Plv ) 

( c.P R b) ( fP L v) 


(l,3)o 

\ 

(gqQLT^qL + gd L T'1 ,I t L )Wl 

Os L 

Ot 

(■ cP L b ) ( fP L v) 
(ca^P L b) {tcf^Plv) 


HR 2 ) 1/2 

(■ X d qLd R (j> ± X u q L u R iT 2 (jA ± A dLe R <t>) 


o[ 


v L 


{r'r^PLb) {c^Plv) 


Ov L 


( 3 , 3)2/3 


°V R 

(f'y fJ ,P R b) (c'y^Plv) < —> 

0’s R 

(■ fP R b ) (cP L u) < —> 

0's L 

{fP L b) (cP L u) < —> 

0' T 

(t(t mi/ Pl6) (cav, v P L v) «■—> 

n" 

u v L 

(f7 M Pnc c ) (bF'fP L u) «—► 

°V R 

(f^P R c c )(b c rPLP) <—*• 

°S R 

(■ fP R c c ) (b c P L v) <—*■ 

0's L 

(fP L c c ) ( b c P L v ) ► 


(fa^P L c c ) (Vv^Plv) <—*• 


-20 Sr 

-- 2 Ov r 


y( 3 ,1)2/3 


-\Os L - | O t 
-6 Os L + \O t 

(3, 2)7/6 

~°v R 


-2 Os R 

(3, 2)5/3 


(3,3) 1/s 

\ 

-\Os L ± | O t 

j (3,1)1/3 

-6 Os L - \O t 



A q L r^ L U^ 

(A qL-y^L + A rfi ? 7 M e J? )?7' i 

(A u R l L + XqLir 2 e R )R 


(Xd^lL + X^e^V^ 
Xq^iT2r£ L S 

{Xq^ir 2 £L + X u R e R )S 


TABLE II. All possible four-fermion operators that can contribute to B — > D^tv. Operators for which no quantum numbers 
are given can only arise from dimension-8 operators in a gauge invariant completion. For other operators the interaction terms 
which are subsequently integrated out are given. For the T operators we use the conventional definition of = ifq'', 7"]/2. 


single operator, while others can generate two simulta¬ 
neously. Anticipating the large Wilson coefficients neces¬ 
sary to fit the observed R(D^) ratios, we focus on oper¬ 
ators which can arise from dimension-6 gauge-invariant 
terms. Operators which can only come from SM gauge 
invariant dimension-8 terms or cannot be generated by 
integrating out a low-spin mediator will be omitted. Such 
contributions would be suppressed by additional powers 
of v/A , or could only arise from strongly coupled NP. 

We calculate the contributions of all operators in the 
heavy quark limit [35]. Subleading corrections Our 
method follows that of Ref. [36], and we rederived and 
confirmed those results. (A missing factor of (1 — m 2 /g 2 ) 
has to be inserted in Eq. (10) of Ref. [Mj.) We use the 
most precise single measurement of the Ha 1 form fac¬ 
tor [8], which equals the Isgur-Wise function in the heavy 
quark limit. 

Higher order corrections are neglected, except for 
the following two effects that are known to be signifi¬ 
cant. For scalar operators a numerically sizable term, 
(tob+TO£)»)/(to&-(- m c ) ~ 1.4, arises from (D* |c 75 &|.B) = 
—q IJ ’(D*\cj IJi "/ 5 b\B)/(m b + m c ). We also include the 
leading-log scale dependence of scalar (here Cs is either 
Cs L or Cs R ) and tensor currents m in fits to models 
where they appear simultaneously, 


Cs(m b ) = 
C T {m b ) = 


a s {ni t ) \ 

a s {m b )J 

a s (rn t ) \ 

a s (m b )J 


-12/23 


4/23 


s(M) 

,(m t ) 


- 12/21 


Cs(M). 


4/21 


s{m t )J 


C t (M). (9) 


For numerical calculations we use a reference scale M = 
750 GeV. The sensitivity to this choice is small, as most 
of the running occurs at low scales, between m b and rn t . 

To test the robustness of our results to O (Aqcd /fn c ,b) 
corrections, we varied the slope parameter of the Isgur- 
Wise function, p 2 , by ±0.2 (motivated by Ref. [33]), and 
found less than ler change in the results. We leave consid¬ 
eration of 0(AQCT>/m Ctb ) corrections for the new physics, 
of the sort carried out for 2HDMs in [39!j [40], for future 
work. 

Figure[l]shows the results of x 2 fits to R(D) and R{D*) 
for each of the four-fermion operators in Table [II] individ¬ 
ually. Here and below, our \ 2 includes experimental and 
SM theory uncertainties, but does not include theory un¬ 
certainties on NP, which are subdominant. Throughout, 
we assume that no new large sources of CP violation 
are present, i.e., we assume that the the phases of the 
NP operators are aligned with the phase of the SM vec¬ 
tor operator. A contribution from the Ot operator or a 
modification of the SM contribution proportional to Oy L 
(or any of its equivalents under Fierz identities) provide 
good fits to the data. The Og L operator can also fit 
the total rates, but it leads to q 2 spectra incompatible 
with observations. The operator Og also gives a good 
fit, which is not apparent from only considering the un¬ 
primed operators. (Note that the measurements of R{D) 
and R(D*) depend on the operator coefficients, because 
the decay distributions are modified by the new physics 
contribution, affecting the experimental efficiencies and 
the measured rates nun]. This effect cannot be included 
in our fits, providing another reason to take the y 2 values 
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FIG. 1. Goodness-of-fit for the coefficients of individual operators from the measured R{D) and R(D*) ratios. Besides the 
fits to the unprimed operators in Table [II] (left), we also show fits to primed operators not related by simple rescalings (right). 
Faded regions for Cs L indicate good fits to the observed rates excluded by the measurement of the q 2 spectrum [2]. Note that 
the x 2 includes experimental and SM theory uncertainties, but not theory uncertainties on NP. 
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FIG. 2. Goodness-of-fit for coefficients of operators which can be generated from dimension-6 operators with fermion bilinears 
having the same SM quantum numbers. The plots show 1-, 2-, and 3cr allowed regions. Approximate regions of parameter 
space excluded by the measurement of the q 2 spectrum [2] are presented as faded regions, as in Fig. |T] 


as rough guides only.) 

As noted earlier, certain mediators can generate two 
contributing operators simultaneously. Fig. [2] shows the 
three such two-dimensional y 2 fits. While any two rates 
can be explained by fitting two operator coefficients, the 
existence of a solution consistent with all other con¬ 
straints with a given flavor structure is nontrivial and 
is the topic of the following section. A summary of all 
coefficients of best fit points with y 2 , ; „ < 5 and accept¬ 


able q 2 spectra is provided in Table III 


Besides the branching ratios, additional model discrim¬ 
ination comes from the q 2 spectra (especially in B —» 
DtD ), which are consistent with SM expectations J21 Sl¬ 
it is not possible to do a combined fit with publicly avail¬ 
able data, because correlations among different q 2 bins 
are unavailable. We follow Ref. j2] in eliminating cer¬ 
tain models by comparing their predicted q 2 spectra with 
the measurement. It was observed that two of the four 
solutions in the Cs R -Cs L plane (Fig. [2| left plot) are 


Coefficient (s) 

Best fit value(s) (A = 1 TeV) 

C Vl 

0.18 ±0.04, 

-2.88 ±0.04 

Ct 

0.52 ± 0.02, 

-0.07 ±0.02 

c'L 

-0.46 ±0.09 


(■ Cr , C L ) 

(1.25,-1.02), 

(-2.84,3.08) 

(Cy R , Cy L ) 

(-0.01,0.18), 

(0.01,-2.88) 

(CS*. CgJ 

(0.35,-0.03), 

(0.96,2.41), 


(-5.74,0.03), 

(-6.34,-2.39) 


TABLE III. Best-fit operator coefficients with acceptable 
q 2 spectra and Xmin < 5- F° r the ID fits in Fig. [l] we in¬ 
clude the Ay 2 < 1 ranges (upper part), and show the central 
values of the 2D fits in Fig. [5] (lower part). 


excluded [2], as indicated by the faded regions. In the 
Cy R -Cy L plane (middle plot), we find the measured q 2 
spectra exclude regions that provide good fits to the total 
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FIG. 3. Normalized dF(B —> Drv)/Aq 2 distributions. The 
histogram shows the BaBar data [2]. The red, blue, and green 
curves correspond to models with Wilson coefficients in the 


(a) top, (b) middle, and (c) bottom lines in Table IV 


couplings are excluded by other rare decays and CKM 
unitarity constraints [47]. We therefore consider NP that 
is charged under a flavor symmetry, focusing on the pos¬ 
sibility of Minimal Flavor Violation (MFV) [4SHSD]. In 
the quark sector, the MFV framework assumes that the 
breaking of the U(3)q x U( 3) u x U(3)d flavor symmetry 
has a single source in both the SM and beyond. It is 
parametrized by the Yukawas, acting as spurions of the 
symmetry breaking, transforming in the (3,3,1) repre¬ 
sentation for Y u and the (3,1,3) for Yd. Focusing pri¬ 
marily on MFV in the quark sector, for viable models 
we attempt to extend the MFV scenario to the lepton 
sector as well. We find a unique model that is consis¬ 
tent with MFV both in the quark and the lepton sector, 
and comment on alternative approaches involving hori¬ 
zontal symmetries. We begin by showing in Sec. |I TTK\ 
that uncolored mediators are disfavored by other con¬ 
straints. Then, in Sec. |IIIB[ we identify viable MFV 
models with leptoquark mediators. 


Label 

Coefficients (A 

= 1 TeV) 

Comment 

(a) 

C' Vr = 1.10 

C' Vl = 0.24 

disfavored 

(b) 

C' Vr = -0.01 

C' Vl = 0.18 

allowed 

(c) 

C'L = 0-96 

CL = 2-41 

allowed 


TABLE IV. Operator coefficients for the q 2 spectra in Fig. [ 3 ] 

rates for values of \C' Vr | >0.5. In the Cg R -Cg L plane 
(right plot) all fits consistent with the total rates are also 
consistent with the q 2 spectra. Figure [3] illustrates this 
by showing the measured dr (13 -4 Drv)/Aq 2 spectrum 
together with one disfavored and two viable models, cor¬ 
responding to the entries in Table m 


III. MODELS 

Having identified dimension-6 operators that can gen¬ 
erate the observed B —> D^tD rate, we now consider 
possible UV completions. The TeV-scale required for the 
four-fermion operators points to a tree-level NP contri¬ 
bution. Concerning the flavor structure, one possibility 
is that NP is aligned with the SM Yukawa matrices and 
in the fermion mass basis only gives rise to the be Dt four- 
fermion interaction. In such a scenario the measurement 
of R(D^) would have little interplay with other flavor 
data; however, it is hard to imagine a UV completion 
with such precise alignment with the SM flavor struc¬ 
ture. On the other hand, if operators are generated with 
general fermion content, Pijki diiij i\ii, where i, j. k, l are 
generation indices, then some Pijki coefficients must sat¬ 
isfy strong experimental constraints. Most of the com¬ 
pletions we consider below have also been studied, but 
without considerations to their possible flavor structure; 
see, e.g., Refs. EHiili 

Independent of the UV completion, flavor-anarchic 


A. Uncolored mediators 

1. Higgs-like scalars 

We first explain why no new color-neutral scalar, 
such as those in the nonstandard 2HDMs proposed in 
Refs. [5TH551 . is compatible with an MFV structure. 
The simplest case, a flavor singlet scalar, would be con¬ 
strained to have couplings proportional to the SM Higgs. 
As shown in Fig. [2] comparable values of Cs L and Cs R 
are needed to explain the current B -4 D^tv data. Cs L 
is proportional to y c , which implies that B —i D^tv 
cannot be fit, keeping the charged Higgs heavier than 
collider limits, consistently with perturbativity of y t . 

An alternative is to charge the scalar under the quark 
flavor symmetries. In this case, the scalar needs to be 
in the representation of some combination of Yukawa 
spurions, since it must couple to leptons as well. To 
generate both scalar couplings simultaneously, an un¬ 
suppressed 0{ 1) coupling to one chiral combination of 
first-generation quarks is unavoidable. Integrating out 
the scalar generates dangerous four-quark and two-quark 
two-lepton operators. Four-quark operators of this type 
are strongly constrained, and such a scalar would also 
appear in t + t~ resonance searches [54, '55]. 

2. W'-like vector triplets 

A rescaling of the SM operator, Oy ,, provides a good 
fit to the data. A simple possibility is then the presence 
of a W', a new vector that couples similarly to the SM 
W boson. To avoid explicitly breaking the SM gauge 
symmetry, a full SU(2)l triplet is needed. The simplest 
choice, a flavor singlet IV', is tightly constrained by the 
LHC [55H55] , with electroweak-strength coupling of a W' 
to first generation quarks excluded up to mw 1 ~ 1.8 TeV. 
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This is in conflict with the 0.2 ~ g 2 \V c b\ (1 TeV/mw') 2 
coupling needed to explain the B —y D^tv data. 

As with the scalars, it is possible to consider a W' 
that transforms under a nontrivial representation of the 
[t/(3)] 3 quark flavor symmetry group, in an attempt to 
suppress the couplings to light quark generations. How¬ 
ever, the W' needs to simultaneously couple to leptons, 
fixing its flavor assignment to be the conjugate of some 
combination of Yukawa spurions. This means that the 
same flavor structure allowing the W' to couple to lep¬ 
tons can always couple to a current of left-handed quarks 
contracted as a flavor singlet, allowing for universal cou¬ 
plings to all quark generations. 

Suppose that the flavor singlet universal coupling to 
quarks is accidentally small. The minimal charge assign¬ 
ments are ( 3 , 3 , 1 ) and ( 3 , 1 , 3 ). In the ( 3 , 3 , 1 ) case, 
(semi)leptonic b —y c transitions are suppressed by y c , 
and therefore the central values of R(D^) cannot be 
generated with perturbative couplings. In the ( 3 , 1 , 3 ) 
case, the R(D ( '*' ) ) data can be accommodated by the in¬ 
teraction SC = gqjY^T • WY‘C‘q\ and an equivalent 
coupling to the leptons with the flavor indices of Y d and 
W M contracted into each other, given yb = 0(1). LHC di¬ 
rect production bounds are then evaded. However, tree- 
level FCNCs are mediated by the neutral component, 
W°, and fitting the B — > tv rates implies too large 

of a contribution to B-B mixing. For an approach us¬ 
ing a combination of dynamical singlet suppression and 
reduced symmetries to evade these issues, see Ref. f59) . 

B. Theories of MFV leptoquarks 


Unless otherwise specified, all comments below apply for 
both cases. Since these leptoquarks need to mediate 
processes involving different generations of quarks, ad¬ 
ditional experimental constraints are important. We will 
see that the dominant bounds come from the tension be¬ 
tween allowing for a large enough effect to fit the R(D^O) 
data, while being consistent with precision electroweak 
constraints and contributions to FCNC decays, such as 
B —y X s vv. 

Most choices of representation under the flavor groups 
can be immediately discarded, since they either only 
generate Yukawa-suppressed operators for the required 
b — y ctv decays, or have unsuppressed couplings to first- 
generation quarks. The latter case would lead to new 
contributions to t + t~ production at the LHC, which is 
severely constrained by Z' searches [531 155) . We estimate 
that couplings larger than A/m ~ 0.25/TeV are excluded 
due to production of t + t~ through t-channel leptoquark 
exchange. Below, we only consider models that satisfy 
both of these constraints. 


1. Scalar leptoquarks — r alignment 

Consider a scalar leptoquark, S , with U(3)q x U(3) u x 
U(3) d charge ( 1 , 1 , 3 ). The interaction terms with the 
smallest number of Yukawa insertions for the electroweak 
singlet case take the form 

SC = S(AY d f <f L ir 2 4 + A YjY u u c R e R ) (12) 

= Si{*y di V*i u c L] t l - A y di d c Li v L + A y^y^V-i u c Rj r R ). 


The second class of potential tree-level mediators are 
leptoquarks. With no way to contract a single quark 
representation with Yukawa spurions to form a flavor sin¬ 
glet, leptoquarks in the MFV framework must carry some 
charges under the quark flavor symmetries for flavor in¬ 
dices to be contracted in interaction terms. Additionally, 
for representations of the flavor group to be coupled to¬ 
gether in a way that gives non-zero couplings, the medi¬ 
ators must come in fundamental representations of some 
of the [U(3)] 3 quark flavor group. Any larger represen¬ 
tation will always either give vanishing coupling when 
contracted with Yukawa spurions in order to couple to 
quarks or have the spurions entirely contracted into the 
leptoquark so as to not generate any new couplings. For 
reviews of leptoquarks, see, e.g., Refs. J60jl6h|, and for 
discussion of MFV leptoquarks, see, e.g., Ref. [621 . 

The six simplest possible flavor assignments under 
U(3)q x U(3) u x U(3) d , for scalar leptoquarks, are 

S'- ( 3 , 1 , 1 ), S- ( 1 , 3 , 1 ), ( 1 , 1 , 3 ), (10) 

while vector leptoquarks can have charge 

Un ~ ( 3 , 1 , 1 ), ( 1 , 3 , 1 ), ~ ( 1 , 1 , 3 ). ( 11 ) 

For each of these six charge assignments, the leptoquarks 
could be either electroweak SU(2)l singlets or triplets. 


Couplings to light quarks are suppressed by quark masses 
and/or small CKM angles, so LHC bounds involving pro¬ 
duction off valence quarks are evaded. Meanwhile, b —y c 
operators receive contributions of the form 


C 


Sr 


Vi 


cb \ 2 2 


A 2 


A Yb 


c 


= isr^ycy'i ■ ( 13 ) 


v: 


l s 3 


m 


S 3 


This contribution can be sizable if yb = 0(1), which can 
follow from a type II 2HDM at large tan/?. Due to the 
large suppression of C Rl from y c , we use from Fig. [2j 
C'l » 0 and C Rr / A 2 = (0.35 ± 0.1)/(1 TeV) 2 , corr^ 
sponding to A yb/ms 3 « (2.95 ± 0.6)/TeV, with A arbi¬ 
trary as long as it remains perturbative. The analogous 
term for the electroweak triplet scalar only generates the 
Cg R term, but with the opposite sign, due to the pres¬ 
ence of an extra r 3 factor in the interaction Lagrangian 
for the relevant component. We do not study the triplet 
further, since fitting the R(D < '*' ) ) data would require a 
large imaginary part for A. (In this paper we neglect 
the possibility of complex couplings; otherwise CP vio¬ 
lation constraints need to be studied, which go beyond 
the scope of this paper and are unlikely to result in qual¬ 
itatively different allowed scenarios.) Similar operators 
contribute to B~ —y tv decay, with V c b replaced by V u b . 
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The contribution of these operators to the measured rate 
is consistent with current data. 


Coupling exclusively to third generation leptons can 
be enforced by imposing a horizontal flavor symmetry, 
U(1) T , under which the leptoquark and third-generation 
lepton fields are equally and oppositely charged. Neu¬ 
trino masses and mixings require additional spurions that 
break the lepton flavor symmetry, but these effects can 
be small enough to ignore in the present discussion. For 
an example where such symmetries generate the entire 
structure of b —► ctv transitions, see Ref. [53], and 
Refs. [64] [65] for discussions of lepton symmetries in 
b — > c transitions in relation to other anomalies. 


The presence of a large S$t c T coupling also generates 
a deviation from the SM prediction for the Zt + t~ cou¬ 
pling, which can affect both the total rate and asymmetry 


of Z 


r^'r 


decays. For scalar leptoquarks, such cor¬ 


rections are finite and calculable |66| . In the case where 
only the couplings proportional to A are present, these 
provide a bound on A yb/ms 3 comparable to the value 
needed to fit the R(D^) data. However, interference 
between the A and A couplings can be large, since the 
latter is not suppressed by y c (unlike in the decay of the 
6-quark). If |A/A| > 0.45, interference among A and A 
contributions can render the model unobservable to cur¬ 
rent electroweak precision tests. 

Being color triplets, leptoquarks can be pair-produced 
directly at hadron colliders from gg initial states, inde¬ 
pendent of their couplings to fermions. The leptoquark 
that contributes to B — > D'*>tv, S3, decays almost en¬ 
tirely to either tr or bv final states, with relative branch¬ 
ing fraction determined by the A and A couplings, but 
with a branching fraction to tr not less than « 0.5. 
Such searches have been carried out in both channels 
by CMS [57]. Assuming equal branching fractions for a 
conservative bound, the mass of the S 3 leptoquark is re¬ 
quired to be ms 3 > 560 GeV. Limits from decay purely 
to the bv final state can also be derived by reinterpret¬ 
ing the bounds from sbottom searches [55] [55], providing 
comparable limits in areas of parameter space where they 
are applicable. 

While the leading MFV interactions do not mediate a 
contribution to b —> svv, the SM rate is sufficiently small 
that we also need to consider effects subleading in the 
MFV expansion. The interactions in Eq. (12) refer to 
the mass eigenstates only if the leptoquarks’ couplings 
to different down-type quarks do not mix. If such mixing 
occurs, it could induce large b —> svv transition. In a 
spurion analysis, we should consider any number of spu- 
rion insertions. Since we need yb to be 0(1), we should 
include spurions of Yd at all orders 123. and the above 
conclusions should be checked after allowing such inser¬ 
tions. Thus, the leptoquark mass matrix takes the form 
(again in the quark mass basis) 


m 1 2 s =ml(I + bY^Yd + • • •) 
=* rn 2 * s . = ml (1 + by}. + •••)• 


( 14 ) 


This makes it clear that the leptoquarks do not mix to all 
orders in the MFV expansion, while the third generation 
leptoquark can be split in mass from the first two by a 
sizable amount. 

Although the leptoquark mass matrix does not induce 
large b —> svv transitions, strong constraints arise from 
interactions at next order in the MFV expansion^] At 
leading order, Eq. (12) couples down-type quarks to neu¬ 
trinos, but as explained above, does not allow them to 
mix. The next term, with three Yukawa spurion inser¬ 
tions, 


6 £' = A'SYjY u Yjg c L iT 2 t L , (15) 

generates interactions of the form 

S£' = S t A' y di V; x yl. (u c L jT L - V jk d c Lk v L ) . (16) 

This interaction generates b —¥ svv decay, to which the 
leading contribution is 

Obsuv = --V t * b V ts (bL'fSL VL W ) , (17) 

with no parametric suppression relative to b —► c. Such 
transitions are most strongly constrained at present 

by |Z5 

B(B + -> K+vv) < 1.6 x 1CT 5 . (18) 


Comparing the SM prediction with this bound d Eg 
gives X'Xy 2 /(2mg 3 ) < 0.26/(1 TeV) 2 . If Cg R in Eq. (13) 
is to explain the values of R(D^), this can be rewritten 
as A'/A < 0.1. Thus, a moderate suppression of the sub¬ 
leading coefficient of the MFV expansion is necessary to 
avoid the b —> svv constraint. 


2. Scalar leptoquarks — lepton MFV 

We can extend the MFV framework to the lepton sec¬ 
tor, and give the leptoquark an appropriate charge un¬ 
der the U(3)l x U( 3) e symmetry broken by the lepton 
Yukawa couplings. Assigning the leptoquark to either of 
the lepton flavor representations ( 3 , 1 ) or ( 1 , 3 ) ensures 
that either the terms with A or A, respectively, will cou¬ 
ple to all lepton generations uniformly. In the first case, 
the leptoquark contributes as much to B —> D^lv as 
to final states with a r, up to effects of subleading spu¬ 
rion insertions, and R(D < '*' > ) cannot deviate significantly 
from the SM prediction. In the second case, all couplings 


1 Due to the large value of yt, terms to all orders in (Y u Yy ) n should 

also be included in the expansion analogously to the (Y^YJ)" 

terms in Eq. fH- Up to terms suppressed by m c /mt, however, 

this only leads to a rescaling of the coefficient in the first sub¬ 

leading term. 








generating corrections to B —> D^*Hv are suppressed by 
Yukawas of either lighter quarks or leptons. 

For this second case, the leading interaction terms take 
the form 

SC = Su(XY i 2jjY em i q c L jiT 2 (- Lm + A Y d jiY U j k u Rk e Rl ) 

= Su(Xy di Ve.Vji u c Lj r Ll - Xy di y ei d c Li u Ll (19) 

+ A yd i yu j v 3 * i u c Rj T Rl ). 

The leptoquarks coupling to lighter lepton generations 
introduce additional constraints. Of those involving the 
A coupling, the leading constraint comes from consider¬ 
ing D° —» y + y,~ decays, which receive a tree-level con¬ 
tribution from t-channel leptoquark exchange. However, 
the contribution to the branching ratio is suppressed by 
0 (yf l ) (in addition to the already-present helicity sup¬ 
pression), making it two orders of magnitude below cur¬ 
rent bounds m when consistent with the fit to R(D^). 
A contribution to B s —> y + y~ is generated by a mixed 
W-S 32 box diagram, but the resulting correction to the 
SM rate is at the percent level. A more severe constraint 
on A exists from corrections to measured at LEP, 
due to the unsuppressed coupling of the St R £ R vertex, 
providing a limit A/A < 0.8. 

The leptoquarks’ coupling to light leptons mean that 
collider searches for first- and second-generation lepto¬ 
quarks become relevant. These bounds are less constrain¬ 
ing then the third-generation searches in Sec. |III B 1[ 
however, since they involve different leptoquark compo¬ 
nents than that contributing to !?(£)(*)). Here, the mass 
matrix takes the form 

( m ! )ijlm = ml(I + bY^ ki Y dk j + • • • + dY* kl Y ekm + • • •) 
=ml{l + by 2 d . +dy 2 ei +•••), (20) 

so the ‘‘5.33 component can be lighter than the others by an 
0(1) amount without a sizable tuning. The fact that the 
bound on doubly-produced first- and second-generation 
leptoquarks are approximately a factor of 2 more strin¬ 
gent than those for the third-generation ones [73175] does 
not then lead to any unavoidable additional constraints. 
Furthermore, the additional leptoquarks do not lead to 
more stringent constraints from b —> svv decays, since 
all couplings to neutrinos other than v T are Yukawa sup¬ 
pressed. 


3. Vector leptoquarks 

Consider a vector leptoquark, W*, with U(3 )q x 
U(3) u x U(3) d charge ( 1 , 1 , 3 ). If the leptoquark is an 
electroweak singlet, the lowest order MFV contribution 
in the quark mass basis is 

SC = (A q L Y d jfj,e L + X dn^^en) U^ (21) 

= (Xy di V ji u L j"/ ll is L + Xy^dLij^TL + Xdmj^Tjt) U?, 


while for the electroweak triplet, it is 


6C = (Xq L Y d f^e L )U» 

= (a y di Vji ULj'y^L + A y di dLi7 M T L ) (22) 

+ ^ + ydiVji ULj'y^TL . 


(Here the signs on U do not indicate charges, but rather 
which component of the triplet we refer to.) The contri¬ 
bution to b —> c decays comes from 

%-%**■ (23) 

As in the prior case, this leptoquark model requires 
yb = 0(1) to be able to explain the R(D^) data with 
perturbative couplings. In this case, it is the electroweak 
singlet that generates an operator with the wrong sign 
to be fit with real values of A. Using Table [TT| and Fig. [I] 
we see that a value of A yb/mu 3 ~ (2.2 ± 0.4)/TeV gives 
a good fit. 

If we consider giving leptonic flavor charges to the lep¬ 
toquarks, as in the scalar case, then the ( 3 , 1 ) representa¬ 
tion prevents a deviation of R(D^) from being generated 
as before. The (1,3) representation does not, but in this 
case, a tree-level contribution to B s — > [771 [ZB] is 

generated at leading order in the MFV expansion. We 
estimate this to give a constraint of Xyb/mu 3 < 1.0/TeV, 
in significant tension with the value required to fit the 
B decay data. However, coupling to purely r can be en¬ 
forced via a horizontal symmetry, as in Sec. |III B 1[ and 
this is what we assume for the rest of the discussion. 

As in the scalar case, we need to consider the pos¬ 
sibility of generating sizable contributions to b —svv 
from additional spurion insertions. The leptoquark mass 
matrix remains diagonal. As in the scalar case, for the 
electroweak triplet vector leptoquark, the next contribu¬ 
tion to the MFV expansion of the interaction term, with 
three Yukawa spurion insertions, 

SC = X' q L Y u YjY d f^e L UC (24) 


leads to an operator of the form 


VIV 


Obsuv = tb 2 tS VtUb A 7 A (bL^SL U7/iU) • (25) 


m 


u 3 


In this case, the bound is X'Xy^/m^ < 0.13/(1 TeV) 2 . If 
the operators in Eq. (23) are to explain the R(D W) data, 


this can be rewritten as A'/A < 0.05. The electroweak 
triplet case thus requires a suppression of the subleading 
coefficient of the MFV expansion to avoid violating the 
b —> svv bound. There is no such constraint for the elec¬ 
troweak singlet vector leptoquark, however, that model 
requires complex couplings to fit the !?(£>(*)) data. 

There are bounds on the leptoquarks produced directly 
through their gauge couplings via gg —> UU pair produc¬ 
tion. The leptoquarks can then decay to br or tv pairs. 
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mu 3 [GeV] 


FIG. 4. Limits on mu 3 from direct pp —> U 3 U 3 production. 


With current constraints [5^, and again fixing X 2 /mfj 3 to 
explain the observed excess, this implies A ' 2 < 0.1, con¬ 
sistent with reasonable values of A without any greater 
tuning than already necessary. 

In the MFV framework, we thus find that two models 
are compatible with the data, once mild tuning is allowed 
for subleading MFV couplings to evade constraints from 
b —> suu: a scalar leptoquark in the (3, 1 ) 1/3 and a vec¬ 
tor leptoquark in the ( 3 , 3) 2 /3 representation of the SM 
gauge groups and in the (anti)fundamental of the U(3)d 
flavor symmetry group. Contrary to assumptions made 
elsewhere, although highly constrained, MFV physics can 
explain the R(D^) data. 

IV. SUMMARY AND FUTURE SIGNALS 


Direct leptoquark searches exist for the first decay mode, 
while bounds on the second mode can be derived by re¬ 
casting stop searches. Vector leptoquark production at 
hadron colliders is complicated by the fact that in the ab¬ 
sence of a UV completion, a marginal coupling not fixed 
by the minimal coupling prescription is present. In addi¬ 
tion to the interactions given by Eq. (21), the rest of the 


Lagrangian for the vector leptoquark is 


Cu = -\ul v U» v +mlUlU»-ig s K{Ult a U v )G^, (26) 


up to higher-order operators suppressed by the cutoff. 
Here U^ — D^U^, while an additional gauge- 

invariant “dipole term” beyond the minimal coupling pre¬ 
scription is also present. If the vector leptoquark is a 
massive gauge boson from a spontaneously broken sym¬ 
metry, then k = 1 , but in principle k can be arbitrary. 
Setting k = 0 gives the minimal production cross section. 
In all cases, a bound is produced directly on leptoquark 
masses, since the leptoquarks are dominantly produced 
due to their gauge couplings, and not their couplings to 
fermions. These are shown in Fig. [4] Since the lepto¬ 
quark direct production rate is higher than that of stops, 
data up to the leptoquark exclusion limit is not pub¬ 
licly available (although stop searches could be recast for 
higher masses, which we leave for future work). In the 
figure we present a conservative extrapolation. 

Additional bounds come from meson mixing and pre- 


were computed in Ref. m ■ The interpretation of these is 
complicated by the fact that massive vector leptoquarks 
are not in themselves UV-complete and yield logarithmic 
divergences at loop level. In addition to b —> suD, the 
same interactions that lead to Eq. (25), also give rise to 


a quadratically divergent contribution to Bd yS mixing at 
one-loop level. As a conservative estimate, the contribu¬ 
tion of the finite part is then 


OAB=2 = M^ Kv ^ CbLraLf . 


l u 3 


(27) 


We studied possible explanations of the observed en¬ 
hancements of the B —>■ D^tv rates compared to the 
SM predictions. We identified which higher dimension 
operators provide good fits to the data and determined 
their required coefficients. Since a sizable modification of 
a tree-level process in the SM is required, the enhance¬ 
ments of B —> D^tv point to a light, tree-level, me¬ 
diator. We found viable models, consistent with MFV, 
where the mediator is a leptoquark. While these models 
are consistent with present flavor and LHC data, the low 
mass scale of the mediator, m < 1 TeV, implies a variety 
of promising signals at future experiments. 

With more precise future data from Belle II and LHCb, 
the necessary size of flavor violation may decrease while 
maintaining a high significance for NP. In this paper, we 
focused on NP that can accommodate the central val¬ 
ues of the current measurements, but in reality NP may 
have a heavier mass scale, alleviating possible tensions 
with present constraints. If the magnitude of the devi¬ 
ation from the SM decreases while its experimental sig¬ 
nificance does not, it will become important to control 
theoretical uncertainties as well as possible. It would 
likely be advantageous to consider ratios (81. in which 
the range of q 2 integration is the same in the numerator 
and the denominator, 


R(D ( * } ) = 


KmB-m^,)) 2 £)(*) 

tv) 

<_ dq 2 

f(m B -m DM f dr ^ B £)(*);p) 
ml dq 2 


dq 2 


dq 2 


(28) 


which have smaller theoretical uncertainties than 
R(D^). The reason is that including the 0 < q 2 < m 2 
region in the denominator simply dilutes this ratio and 
the sensitivity to new physics effects, and also because the 
uncertainties of the semileptonic form factors increase at 
smaller q 2 (larger recoil). 

We conclude by enumerating future experimental mea¬ 
surements and new physics search channels that can test 
for models motivated by the measured B —> D^tu rates. 













10 


1. More precise data on B —> D^*Wv, including mea¬ 
surements of differential distributions and their cor¬ 
relations. These would also help discriminate NP 
from the SM, and between different NP models. 

2. D —> irvv could be near the 10 -5 level, which may 
be observable in BES III [82]. 

3. In some scenarios a nonresonant contribution to 
t —> ct + t~ is possible. While current t —> cZ 
searches (83j [84] focus on mi+i- consistent with 
mz, there could also be a contribution that is in¬ 
consistent with the Z lineshape. 

4. A NP contribution to 1 -> brv of a few times 
10~ 3 r t is possible. Noting that the tv do not 
come from a W decay in the NP contribution, a 
high transverse mass tail, ?> mw , can 
be searched for in semileptonic t — > t decays. 

5. In this paper we neglected CP violation. A simple 
experimental test is to constrain the difference be¬ 
tween B + and B~ rates to D^tv] a measurable 
effect would be a clear sign of new physics. 

6 . Models that contribute to B —> D^tv can mod¬ 
ify Z — > t + t~ at the order of present constraints. 
Therefore, future precision Z -pole measurements 
may detect deviations from the SM. 

Finally, we note that there may be viable models in 
which the B —> D^tv rates are given by the SM, but 
B — >• D^lv (l = e,fi) are suppressed by interference 
between NP and the SM [55]. 
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Appendix A: Comment on U(2) 3 scenarios 

We presented several MFV models capable of describ¬ 
ing the data while also coupling to other quark 

generations consistently with all constraints. As detailed 


earlier, extending the MFV framework to the lepton sec¬ 
tor in these scenarios is not always possible, and so we 
have assumed horizontal symmetries governing the struc¬ 
ture of lepton flavor where necessary. However, it is pos¬ 
sible to treat quark and lepton flavor uniformly in all 
scenarios, at the cost of a slightly less predictive flavor 
framework. 

A minimal self-consistent extension of the MFV frame¬ 
work exists in the form of the U( 2) 3 flavor symmetry ap¬ 
proach mm- Here, the UV physics is only assumed 
to preserve a U{2)q x U[2) u x U{2)d symmetry for two 
generations in the quark sector, which is closely aligned 
with the first two generations of the SM. Three spurions 
are now necessary to minimally characterize the quark 
sector flavor symmetry breaking, here denoted as 

A u ~(2,2,l), A d ~ (2,1,2), V ~ (2,1,1), (Al) 
while the quark fields decompose as 

q L ~( 2,1,1), u R ~ (1,2,1), d R ~ (1,1, 2), (A2) 

and the singlet third generation fields are q 3 L , t R: and 
b R . The extension to the lepton sector, with an addi¬ 
tional V(2)\ x U( 2) e symmetry, is straightforward, and 
only a single spurion, A e ~ (2, 2) is needed (ignoring 
m„, which requires additional but much smaller spuri¬ 
ons). Unlike the case of MFV, none of these spurions 
have 0 ( 1) components, so higher-order spurion insertions 
are highly suppressed. 

Coupling predominantly to third generation quarks 
and leptons can now be imposed by demanding that the 
mediator is a singlet under all flavor symmetries, as long 
as currents of light quarks and leptons cannot be written 
as flavor singlets. 

As an illustrative example, we consider a scalar lepto- 
quark. Due to the V spurion, this leptoquark can couple 
to all quark generations. The interaction terms are 

sc = S (Ai <f 3L iT 2 £ 3L + X 2 V q c ir 2 £sL 

+ X 1 i c R T R + X 2 V^A u u c T R ). (A3) 

Coupling to the other lepton generations is impossible, 
since no combination of spurions can absorb their flavor 
charge. The contribution due to terms proportional to A 2 
can be neglected, as it always comes with multiple spu¬ 
rion suppression factors. Then the leading contribution 
by spurion counting to b —► c transitions is given by 

C" VI Cl VI - 

s± = -cbx a X 1 , = X a X lVc . A4 

m s A z nig 

Here A a is in the range (A 1; X 2 ) with the precise value 
set by the coefficients of V in up- and down-type Yukawa 
couplings. In the limit of Ai = A 2 (and Ai = A 2 if we in¬ 
clude higher order terms in our Yukawa diagonalization), 
we recover the contribution of the third generation lepto¬ 
quark of an MFV scenario. This is as we would expect. 
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If we were to add a leptoquark in, say, the ( 1 , 1 , 2 ) repre¬ 
sentation to the model above, setting all the suprion coef¬ 
ficients of terms with i 3 i to A and those with r^/e^ to 
A, we precisely recover the particle content and spurion 
structure of the ( 1 , 1 , 3 ) MFV model. 


and using \V c b\ < 1, 


A 


a. 


\i — Xt\\ 

x b - x t 


(AS) 


where x bt t are free parameters defined in Ref. [8H1 . How¬ 
ever, in this case contributions to b —» svv are generated 
by the same interactions. Similar to Eq. (A5) above, the 
coefficient of this interaction can be written as 


Vcb . A 2 — x b Xi 

o Ai-, 

2 x b - x t 


(A6) 


The analysis of the constraints follows the MFV case. 
The treatment of contributions to b —> svv is different, 
however. The term A Q can be written, ignoring phases 


so one can reduce the rate of the interaction by tuning 
Xb ~ A 2 /A 1 without affecting the b —» ctu rate. A slight 
tuning is then again necessary to avoid this constraint, 
although arising in a different set of couplings. 
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